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From a simple model for the driven motion of a planar interface under the influence of a diffusion 
field we derive a damped nonlinear oscillator equation for the interface position. Inside an unstable 
regime, where the damping term is negative, we find limit-cycle solutions, describing an oscillatory 
propagation of the interface. In case of a growing solidification front this offers a transparent scenario 
for the formation of solute bands in binary alloys, and, taking into account the Mullins-Sekerka 
instability, of banded structures. 
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The interaction of propagating extended defects with 
a diffusion field frequently leads to oscillations or jerky 
motions of the defects. A prime example of such an effect 
is the oscillation of a solidification front, induced by the 
diffusion of the solute component in a dilute binary alloy, 
which is growing in the setup of directional solidification. 
In a large number of metallic materials this leads to the 
formation of banded structures pQ , reflecting a periodic 
array of layers with high and low solute concentrations 
where the former ones show a dendritic microstructure. 
The appearance of similar banding effects has recently 
been discussed in rapid solidification of colloids. 

Layered structures are also generated by the oscillatory 
nucleation of a solid phase under the action of a diffusion 
field [3] . A related phenomenon is the oscillatory zoning, 
observed in solid solutions [3] and in natural minerals [5] . 
Another notable scenario is that of diffusion-controlled 
jerky motions of a driven grain boundary 6 . A similar 
behavior of dislocations in metallic alloys leads to the 
Portevin-Le Chatelier effect [7j , denoting the appearance 
of jerky plastic deformations. We, finally, mention the 
oscillatory motion of a crack tip, which is coated by the 
nucleus of a new phase [8] , replacing the attached cloud 
of a diffusion field. 

Theoretical discussions of such effects either are of a 
phenomenological type, like those in Ref. [7], and partly 
in Refs. pQ and [2], or they rely on a Fokker-Planck [6], 
or a diffusion equation with non-equilibrium boundary 
conditions [5] ■ In all approaches the source of oscillatory 
defect motions is identified as an unstable regime where a 
reduction of the driving force leads to an increase of the 
defect velocity. Additional information is provided by 
simulations, based on phase-field models for directional 
solidification [IU] and for nucleation [3] processes. 



In the present Letter we introduce an extremely simple 
but powerful model for the diffusion-induced oscillatory 
motion of a planar interface, using the language adapted 
to the directional solidification of a dilute binary alloy. 
A major advantage of our approach is that it allows a 
transparent and, to a large extend, analytical evaluation. 
This includes a readjustment of the stability analysis by 
Merchant and Davis [TT] who discovered an oscillatory 
instability, similar to that, discussed earlier by Coriell 
and Sekerka [12] . Also included is a clarifying analysis of 
the so far barely understood low-velocity sections of the 
cyclic trajectories, identified by Carrard et al. Q], and 
by Karma and Sarkissian [5]. The limit-cycle behavior, 
describing the oscillations of the interface deep inside the 
unstable regime, is, finally, in remarkable agreement with 
the simulation results by Conti [TO] . 

Due to the restriction to a planar interface, our model 
is a one-dimensional version of the capillary-wave model, 
derived in Ref. 15 from a phase-field model. It is given 
in dimensionless form by the set of equations 

H=lj dz[C(z,t)-U(z-Z(t))] 2 , (1) 

F = F P - m 2 [Z(t) - v P t] 

for the interface position Z(i), and for the excess solute 
concentration C(z, t) relative to its value Cs = 1 in the 
solid phase. The parameter 7 measures the miscibility 
gap AC — Cl — Cs where Cl is the solute concentration 
in the liquid phase, and p measures the mobility of the 
interface. From the equilibrium condition SH/SC = 
it follows that U{z — Z) is the equilibrium-concentration 
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profile at some fixed temperature Tg . This profile reveals 
the smooth solid-liquid transition region of the original 
phase-field model, and is regarded as an input quantity 
of the model ([I]). It effectively comprises non-equilibrium 
effects of sharp-interface descriptions, which are crucial 
for the behavior in the rapid-growth regime, including 
the solute-trapping effect [T4"] . 

The driving force F includes two quantities, appearing 
in the simplest scenario of directional solidification. One 
of them is a constant temperature gradient S, entering 
the parameter m 2 = A^S/Tm where, in physical units, £ 
measures the width of the interface, visible in the profile 
U(£), Tm is the melting temperature of the pure solvent, 
and A is a numerical pre-factor of order one. Adopting 
from Ref. [TU] typical values for S,Tm,£, the parameter 
m 2 is of order 10 -5 . An independent second quantity is 
the velocity vp, applied in pulling the growing crystal in 
opposite direction to the temperature gradient. The local 
temperature Tp at the steady-state position Z(t) — vpt, 
finally, determines the fragment Fp = A(Tg — Tp)/Tm 
of the driving force F. 

The resulting equations for Z(t) and C(z,t) read 



- Z(t) =F P - m 2 \Z(t) - v P t] 
P 



(2) 



7 / dzU'{z - Z(t))[C(z,t) - U(z - Z{t))] , 



(d t -d 2 )C(z,t) = -U"(z-Z(t)) . 

For steady-state boundary conditions C(z = ±00) = 0, 
they have the stationary solutions Z(t) — vpt, and 
C{z, t) = Cp{z — vpt;vp), resulting in the relations 



-vp = Fp + G P {vp)-Gp{0) , 
P 

r+00 

G P (vp) = - 1 dCU'(OCpiC,vp) , 



(3) 



Cp((;v P ) 



C 



dc'u'(c')eMMC'-0] 



We are primarily interested in the late-stage behavior 
of non-stationary solutions Z(t), and, therefore, look for 
a solution C{z, t) of the last equation in Eqs. ([2]), obeying 
the boundary condition C(z, —00) = 0. This leads to the 
expression 



C(z,t) 



t /*+oo 

dt' / dz' d z >G(z- z',t-t') 
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U'{z> Z{t')) , (4) 



involving the Green function 



Q(z,t)=J ^exp(-£; 



-hrt + ikz) 



(5) 



After the substitutions £ = z — Z(t), = z' — Z(t'), and 
expansion of Z(t') around Z(t), we obtain 



g(t-C + z(t)-z(i/) ) t-i/) = 



(6) 



dk 
2^ 



exp [~k 2 {t - t') + ik(( - C) + ik(t - t')v(t)] 



■ exp 



n>2 



where v(t) = Z(t) = v P + h(t). 

If one temporarily applies the scaling transformations 
h —> m~ 2 h,d t —> m 2 d t , one observes that, whereas v(t) 
remains unchanged, a factor m 2 ™~ 2 is attached to the 
contributions oc d™~ 1 v(t). Therefore, with increasing n 
these terms are progressively negligible in Eq. ([6| due to 
the smallness of m 2 . Neglecting all terms of order n > 2, 
we encounter the quasi-stationary approximation, which 
is often used in phenomenological approaches. As we 
shall see, however, a proper understanding of oscillatory 
motions of the interface requires to incorporate the term 
of order n = 2. 

Evaluation of Eqs. Q - ([6]) then leads to the expression 

C(z,t) = Cp((;v) + vl^-{Cp((;v) + Cp((;0)} (7) 
2 ov z v 

with Cp(C;wp) determined by the last line in Eqs. d3b . 
Insertion of this result into the first equation in Eqs. (pi) 
yields the relation 
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F P - m 2 (Z - Z f 



(8) 



1 d 2 1, 



+ G P (v) - Gp(0)+v-j—-[G p (v) + G P (0)} 

with Gp(v) following from the second line in Eqs. 
For v ~ vp the result Q consistently reduces to the first 
equation in Eqs. (|3j). Subtracting the latter from Eq. (|8j), 
we, finally, find for the displacement h(t) = Z(t) — Zp{t) 
the differential equation 

M(h(t)) h{t) + R(h(t)) + m 2 h{t) = (9) 
where we have introduced the mass and friction functions 



Gp{vp + h) + G P (0) 



1 d 2 

Z OV p Vp 

R{h) = -h- Gp(v P + h) + Gp{v P ) 
P 



(10) 



Eq. Q has the typical appearance of a nonlinear damped 
oscillator, and represents one of the central results of the 
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in 
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G P (v) 



v + 2 



(13) 



which determines all terms in the oscillator equation (|9| . 

The resulting numerical solutions for h(t) above and 
below the Cahn threshold vc are shown in Fig. 1 where 
the approach to a limit cycle in the unstable regime is 
clearly visible. For small distances \vp — vc\/vc <C 1 the 
envelopes of these curves can be calculated analytically 
by the Bogoliubov-Mitropolsky method [15] . To leading 
order one finds 



h(t) = a(t)cosijj{t) 



(14) 
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FIG. 1: Numerical solutions h(t) for 7 = 0.01, p = 100, 
m = 0.003, approaching the value h = for vp — 0.522, 
and a limit cycle for vp — 0.520. 



present Letter. We mention that, due to the singular 
dependence of M(h) on vp + h, the differential equation 
Mm is only valid below the crossover line m 2 oc (vp + h) 3 . 



In order to check the stability of the obvious solution 
h(t) = 0, we linearize Eq. ([9| in h(t), which, due to the 
definitions (10) and the first line in Eqs. ([3]), yields 



M(0) h + F'p(v P )h + m 2 h = 



(11) 



In this ordinary oscillator equation the friction coefficient 
can change sign at some critical velocity vc, defined by 
F' p (vc) = 1/p — G' P (v c ) — 0. A similar stability limit 
has been found by Cahn in grain-boundary motion [BJ. 

For quantitative discussions of the behavior of h(t) we 
now adopt the specific model 



L 7(C) = e(-C)expC + e(C)[2-exp(-C)] , (12) 

derived in Ref. [T3] from a double-parabola phase-field 
model. Then, solving the integrals in Eqs. ([3]), one finds 



where ip(t) is a rapidly oscillating phase, and a(t) is an 
amplitude, obeying the differential equation 



da 
~dt. 



-pi a- p 3 a 



(15) 



with pi = ri(vp — vc), and the parameters ri,p 3 fixed 
by the values of 7, p. The solution of Eq. (151 reads 



a(t) = a 



1 , P3 2 
1 + -O0 



exp [2pi t] — — a 2 . 
Pi 



-1/2 



(16) 



which for pi > and p\ < describes the envelopes in 
Fig. 1. The asymptotic value of the limit-cycle amplitude 
shows the critical behavior a(oo) = \/—pi/p3- 

The numerically obtained limit-cycle trajectories of the 
quantities h(i), h(t) , C (Z (t) , t) deep inside the unstable 
regime are displayed in Fig. 2. They also inform on the 
local temperature at the oscillating interface, since this 
is measured by the quantity m 2 h(t). The pronounced 
oscillations of the solute concentration at the interface 
C(Z(t),t) reflects the appearance of solute bands. From 
the curves in Fig. 2, which in part are remarkably close 
to the findings by Conti in Ref. [TO], one concludes that 
the high- and low-concentration layers are connected by 
large-acceleration segments, explaining the sharpness of 
the interfaces between these layers. 

In order to explore the possible formation of dendritic 
ripples, one has to consider perturbations of the form 
h(x, t) = h(q, uj) exp (iq • x + cut) in a three-dimensional 
version of the model ([!]) . In view of the almost stationary 
behavior of h(t) in Fig. 2 at low velocities, we choose, 
as an approximation, Z(t) = vp as a reference velocity. 
Following Ref. [13], we then find the dispersion relation 



OJ 



+ q 2 + m 2 - vp[G P {vp + A) - Gp(v P )} 



X 2 -g 2 
v P + 2A 



[Gp(vp + \) + Gp(X)} 



(17) 
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FIG. 3: Neutral stability curves, enclosing the regions of 
the Cahn (solid line), and of the Mullins-Sekerka instability 
(dashed line). The vertical lines are projections of limit cycles 
at 7 = 0.01 and 7 = 0.02, both at p = 100, m = 0.003 and 
v P = 0.5. 



metallic alloys. The other small cycle in Fig. 3 generates 
layers of precipitation-free periodic solute concentrations. 

The most obvious generalization of our procedure is 
to explore the formation of non-planar layering effects, 
which also is a field for experimental investigations. A 
typical example of such an effect is the oscillatory growth 
of a spherical nucleus, which has been discussed on the 
basis of a phase-field model in Ref . [3] , and which we are 
going to reconsider within our approach. 
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FIG. 2: Trajectories of h(t),h{t),C(Z(t),t) in the unstable 
regime for 7 = 0.02, p = 100, m = 0.003, and v P = 0.5. 



with A = —(vp/2) + (vp/2) 2 + uj + q 2 where the term 
m 2 is the only new element. The wave-number threshold 
q c for the Mullins-Sekerka instability ;16j is determined 
by the relations u>i(q c ) = u)'i(q c ) = 0. By elimination of q c 
from these equations one generates the neutral-stability 
boundary of the instability in form of a function up (7), 
with a parametric dependence on m. 

In Fig. 3 the projection of the limit cycle, belonging 
to Fig. 2, enters the Mullins-Sekerka unstable regime 
at low velocities where the interface develops a dendritic 
microstructure, a typical feature of banded structures in 
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